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The properties of Green's functions and various correlation functions of density and spin operators 
are considered in a homogeneous spin-1 Bose gas in different phases. The dielectric formalism is 
worked out and the partial coincidence of the one-particle and collective spectra is pointed out below 
the temperature of Bose-Einstein condensation. As an application the formalism is used to give two 
approximations for the propagators and the correlation functions and the spectra of excitations 
including shifts and widths due to the thermal cloud. 
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I. INTRODUCTION 

The recent realization of BEG in an optical trap has opened a new area of research ||-|^ . The theoretical investi- 
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md a number of other problems [|18|-p7[ . Gompared to the BEG realized in magnetic traps ||2^,^ , where the 
spin degree of freedom of the particle field is frozen, such systems, since the spinor nature of the particles is preserved, 
have a wider variety of excitations including spin density waves, transverse spin waves or quadrupolar spin waves. 

It is a well established property of scalar Bose-Einstein condensed systems that the one-particle and the density 
correlation function spectra coincide. To treat this problem consistently the dielectric formalism has proved to be 
particularly useful, which has been worked out first for homogeneous systems ||30|~|35t and recently generalized and 
applied to trap systems |^-^. In this paper this formalism is extended to gases with a spinor Bose-Einstein con- 
densate. To make the presentation more transparent only homogeneous systems will be considered. The new feature 
is the appearance of correlation functions including spin fluctuations besides the density autocorrelation function and 
a variety of new one-particle Green's functions. Their perturbation series are analysed simultaneously. After suitable 
rearrangement of the expansions it has been achieved that certain one-particle Green's functions and correlation func- 
tions have common denominators leading to the coincidence of their spectra, though with different spectral weights. 
Examples are the Green's functions corresponding to a spin transfer zero and the density correlation function, fur- 
thermore those Green's functions and correlation functions of the spin operators, which can be characterized by spin 
transfers -1-1 or —1. Exceptions are found in the polar phase (occurring when the interaction in the spin channel 
is repulsive), namely those correlation functions of the spin operators which create zero or ±2 spin transfers do not 
possess condensate induced intermixing with the one-particle Green's functions. The general theory is illustrated in 
the Bogoliubov theory (valid at very low temperatures) and in the random phase approximation (RPA). Within the 
RPA we do not include exchange processes, that already in the case of the scalar Bose-Einstein condensed systems has 
led to involved calculations [ p9| . Their extension to spinor Bose-Einstein gas will be presented in a separate paper. 
As a matter of fact the relative importance of the exchange processes is less in the present case since the direct terms 
are enlarged by summations over spin variables. 

Though the RPA (as a mean field theory) looses its validity near the phase transition it is enlightening to investigate 
it in the transition region. It turns out that, while the transition to the polar phase is continuous when decreasing the 
temperature, it becomes weakly first order for ferromagnetic ordering due to the small (for this coupling attractive) 
spin dependent part of the interaction. It is worth recalling that in case of the scalar Bose-Einstein condensation the 
transition is second order in this RPA-Hartree model, but including exchange contributions makes it to a first order 
one and the effect is then not small, see Ref. and references therein. 

The paper is organized as follows. In the second section, after giving the specifications of the Hamiltonian and 
introducing the canonical transformation to induce the Bose-Einstein condensation, the normal and anomalous one- 
particle Green's functions (matrices in the spin variables) and the different correlation functions of the particle number 



density and spin density operators are defined. Section III is the backbone of the paper. It starts with a summary 
of symmetry properties of the functions introduced in Section |l[ Then their general structure is analyzed within 
the framework of perturbation expansion. We proceed by classifying the processes according to the spin transfers 
involved. Then in the spirit of the dielectric formalism the proper irreducible graphs are separated. The main results 
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are the description of condensate induced intermixing of one-particle and collective modes and the specification of 
the conditions when it occurs. Sections |^ and ^ are devoted to approximate calculations. In Section corrections 
consisting of damping terms and frequency shifts as well are given in RPA to the Bogoliubo y ap proximation discussed 
in Section We present in Section also results for static properties in RPA. Section |Vl| contains the summary 
and further discussion. 



II. BASIC EQUATIONS AND DEFINITIONS 

A. The effective Hamiltonian of the spin-1 Bose gas 

In a gas containing spin-1 particles the wavefunction is a three-component spinor. To be more specific we use the 
basis set of eigenvectors in the spin space, which leads to the wavefunction 
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and the corresponding raising and lowering operators: 

F+ = V2 
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We consider a system of spin-1 particles in a box with periodic boundary conditions and without an external 
potential. The interaction between the particles is an s-wave scattering. Such a rotationally invariant scattering 
interaction can be described by a potential 



y(ri - rs) = (5(ri - r2)[c„ + c,Fi • Fa], 



(4) 



where the constants c„ and Cs is connected to the s-wave scatteriiig lengths (go and ga) in the total hypcrfinc spin 
channel zero and two with c„ = {go + 2g2)/3 and Cs = (52 — 5o)/3 [ |l4[ . 

In the second quantized formalism we introduce annihilation and creation operators ar(k) and aj(k) which destroy 
and create one particle states of plane waves with momentum k and spin projection r. These operators are bosonic 
in our case, so they satisfy the commutation relations: 



[a,(k),at(k')] = J,,,,(5(k-k'), 
[a,(k),a,(k')] = [4(k),at(k')] =0. 

In this formalism the grand canonical Hamiltonian of the system reads as: 

W = ^(ek-A*)4(k)a.(k) + i ^ (ki)at (k2)F;/'a,(k3)a,, (k4 



(5a) 
(5b) 

(6) 



ki+k2 = 

= k3+k4 



where ek = fi^k^/2M stands for the kinetic energy of a particle and /i is the chemical potential of the system and 

V;/ = cJrsSr's' + C,(F),,(F),,,, (7) 

according to (Q). Here and in the following the convention of summing over repeated indices is applied except when 
stated otherwise. 

There are two types of systems depending on the sign of Cg- If it is less than zero then the spins prefer parallel 
alignment which leads to a macroscopic magnetization in the presence of a Bose-Einstein condensate. If Cg is greater 
than zero then the energetically favorable state is when (F) = 0. The former case is called ferromagnetic case and 
later is called polar case [pj. 
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B. Description of the symmetry breaking 



In the Bose-Einstcin condensed phases some field operators have anomalous averages reflecting a broken gauge 
symmetry: (ar(0)) = \/N^Cr and ^aJ(O)) = V^oCr where Nq is the number of particles in the condensate and is 
the normalized spinor of the condensate |l4[ . For the polar case one can take (^r as (0, 1, 0)"^ and for the ferromagnetic 
case can be taken as (1,0,0)^ where the superscript T denotes the operation of transposition. The averaging 
is made over a symmetry breaking ensemble. To consider this symmetry breaking, one can introduce a new set of 
bosonic annihilation and creation operators with a canonical transformation: 



6t(k) =at(k)-(5k,ov^C^• 



(8a) 
(8b) 



The relation between the chemical potential and condensate density can be derived from the requirement that (6r(k)) = 
(6j(k)) = Here and from now on the averages will be made over the grand canonical ensemble with a density matrix 
p = Z^^e*^^ where Z ~ Ti e^^ is the grand canonical partition function. Substituting the canonical transformation 
into (^) one obtains the Hamiltonian in terms of the new operators: 

n = ^(Ck - flo)bl{k)br{k) - fl^MiCbriO) + bl{0)Cr] - M^O + ^ J] 6^ (ki )6t (ks ) V;,'^' 6, (kg)^,' (k4) 
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^bl,io)Cv;:^' CAs' + %l4v;/CsCs' + E(mo - M^)brik). 



(9) 



Here we added and subtracted the term f-'-obl{k)br{k) {^q < 0) to avoid the difficulty of being /i positive in the 
condensed phases, which would lead to a singularity in the unperturbed propagator. 

It should be stressed that with the canonical transformation (^) one defines bare quasiparticle states Y[i bl.{ki) |0). 
The time evolution of these quasiparticle states are determined by the Hamiltonian which contains terms with 
different number of creation and destruction operators leading to the nonconservation of the total number of quasi- 
particles. 



C. Green's functions and correlation functions 



First we define the Green's functions as 

e;|(k,r) =-(T467(k,r)6f (k,0)]), (10) 

with 

The Greek indices introduced here are for distinguishing between the normal and anomalous Green's functions. 
Automatic summation over repeated Greek indices is understood as for the Roman ones. Here and from now on r is 



3 



the imaginary time and is the r ordering operator 1 40 1 . In this symmetry breaking system because of the fact that 
the Hamiltonian Eq. contains terms with two creation or two annihilation operators anomalous Green's functions 
arise (with "f5 — —1). The expression ( |To| ) is the generalization of the well known normal and anomalous Green's 
functions in case of a complex scalar field [[40y30[ . The propagators ([l0| ) are periodic in r with period j3h so their 
Fourier series can be defined as 

g^HKiu^n) = ^^'"dre-"-t;;|(k,r), (12) 

where LOn = 2mT/ f3H is the Bose discrete Matsubara frequency. The spectrum of the one-particle elementary excitations 
can be determined as poles of the analytic continuations of the Green's functions. 

The collective excitations in the system can be described with the correlation functions of the following operators: 



n(k)=^at(kH 
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Here n(k) is the particle density operator .^^(k) is the spin z component density operator, J^±{\i) is the density 
operator of the ±1 spin raising or lowering operator and J^g(k) is the density of the ±2 spin raising or lowering 
operators and at last CTrs(k) is the general density operator from which the other ones can be easily calculated as e.g. 
n(k) = (T+_|_(k) + croo(k) + cr__(k) or J-z(k) — cr++(k) — cr__(k). The different collective excitations can be found as 
poles of the analytical continuations of the corresponding correlation functions (defined for k 7^ 0) , given by 



D„n(k,T) = -(r,[n(k,T)nt(k,0)]), (14a) 

D.^iKr) = -(Tr[T,{k,r)J^l{k,0)]), (14b) 

i?„,(k,T) = -(Tr[n{k,T)J^l{k,0)]y (14c) 

D±±ik,T) = -(^Tr[T±ik,T)Tl{k,0)]), (14d) 

D2±{k,T) = -(Tr[T2{k,T)T2Hk,0)]). (14e) 



Note that n'l'(k) = n(-k), .F|(k) = J^^(-k), J^i(k) = T^{-k) and J^±'^(k) = J^^(-k). A general correlation function 
can also be defined as 



D:.r,,{k,T) = -(^Tr[<Jrs{k,T)as'r'{-KO)] 

from what the above ones can be calculated as: 

A.„(k,T)=^i?-(k,r), 
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All of these correlation functions are also periodic in r with period j3H, so their Fourier transforms can be defined 
analogously to Eq. (|l^) and the appearing Matsubara frequencies are also the same as for the Green's functions. 
We will also apply anomalous correlation functions 

A- (k, r) - - (Tr (k, T)6f (k, 0)] ) , (16a) 

A^^fAKr) ^ -{Tr[b:(k,T)as>r'{^KO)]), (16b) 

which have zero value in the symmetric phase (i.e. for a noncondensed system). 

Though the one-particle Green's functions are the autocorrelation functions of the order parameter field operator, 
we preserve in the following the term correlation function for the other correlation functions introduced above. 



III. GENERAL FORMALISM 



In this section the general formalism is worked out. First symmetry properties of the functions defined in the 
previous section will be discussed. Next perturbation theory will be presented and then we will turn our attention to 
relations valid in all orders of perturbation theory. 



A. Symmetry properties 

Due to the rotational symmetry in the coordinate space the Green's functions ( p^ and all the correlation functions 
( |l^ ) and (|l^) depend only on the modulus of the momentum. 

For the Green's functions the following symmetry properties can be derived. Since (in absence of a magnetic filed) 
the Green's functions (|l^) are real and the cyclic property of trace 

GI^^siK r) = -It, [67 (k, T)6f^" (k, 0)] ) = -{Tr [KiK r)bf (k, 0)] ^ ) 

^-{Tr[biik,r)bfik,0)])=gi;ik,r) 



holds. Furthermore the time displacement symmetry of Hamiltonian systems, the 67(k) = 6,, '''^(— k) relation and that 
after the order of the bosonic operators are irrelevant lead to 

ty;|(k,r) =-(r4&7(k,r)fef (k,0)]) =-(T46f (k,0)fe7(k,T)]) 

= -(T,[bf{k,-T)bUk,0)])=-(Tr[b:'{-K-r)b;''\-KO)]) (18) 

= gLl_^{-k,-T). 

On the basis of the above equalities one finds 

f;;Kk,T) = ty;Kfc,T) = ty|!;(fc,T) = g^_i_^{k,-T) = t;- ^_,(fc,~r). (19) 



For the generalized density correlation functions of Eq. (14i) the following symmetry relations can be derived: 
-i^^r,,(k,T) = (r4a,,(k,r)a,v'(-k,0)]t) = (T,[a^V(-k,0)at,(k, -r)] ' 



(20) 

= {Tr[ar's'ik,T)a,r{~k,0)]) = -Di^' ik,T) 



and 

-D-:.r,,{k,T) - (Tr[ars{k,T)as>r'{-k,0)]) = (Tr[<Js'r'{-k,0)ars{k,T)] 



(21) 

Tr[as'r'{-k,-T)arsik,0)]) ^ - Dl'/ {-k, -t) , 



where we further used that crtj,(k) = asri—k). Combining Eqs. ( pO| ) and ( ^l| ) together one arrives at 

D^^^AKr) = D:.:Ak,r) = D^^' {k,T) = D:'/{k,-T) = I?- ,(fc,-r). (22) 
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For the anomalous correlation functions the two symmetry properties can be derived in the same way. The first 
one cames form the fact that the expectation values here are also real 

A:;(k, r) = ~(Tr K.(k, T)bf (k, 0)] ^) = -(t. [6^(k, tM-K 0)] ) = A^^k, r). (23) 



The second one can derived using the invariance of the trace under cyclic permutations and that the order is irrelevant 
behind an ordering operator: 

A- (-k, -r) = -(r, [<Jrs{-k, -r)6f (-k, 0)] )^-{Tr [&f (-k, r)a.,(-k, 0)] ) = A:r"(k, r). (24) 

Combining Eq. ( |2^ ) and Eq. ( p4| ) together results in (also using that the momentum dependence comes from k = |k|) 
AZ^{k,T) = AZ,{k,T) = A^^f{k,T) = A':r''(k,-T) = A::_^{k,-T). (25) 

B. Perturbation theory 

To calculate the averages of the grand canonical ensemble with the Hamiltonian (|^) one can use the methods of the 
finite temperature many body physics pO| . 

The Hamiltonian of the non interacting system is the first term in Eq. (^): 

Wo = ^(ek-Mo)&]:(k)5,(k), (26) 

k 

which defines free Green's functions: 

^^o^),,(k,z^„) = (27) 

This will be symbolized as a line or a line with an arrow if the Greek indices are specified as seen in Fig. 0. The full 
Green's functions will be symbolized with double lines as shown in Fig. ^ The terms in the Hamiltonian Eq. (|^) 
without any filed operator can be disregarded for our purposes. Among the interaction terms there is one containing 
four field operators and corresponding to a scattering of two quasi-particles which are noncondensed before and after 
the collision as well. There are four interaction terms containing three field operators and six terms containing two 
field operators, corresponding to scattering processes involving both condensate and noncondensate atoms. There 
are interaction terms containing one field operator describing scattering processes involving three condensate and one 
noncondensate atoms. These interaction terms will be graphically represented as shown in Fig. ^. There is one 
remaining term with one field operator and the last term with two operators corresponding to vertices with one and 
two legs as seen in Fig. 0. 
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FIG. 1. The graphical representation of the free propagators 
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FIG. 2. The graphical representation of the full propagators 
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{-h-')NoCv;/C' , {-h-')Nl'\lxlvf/ Cs 



FIG. 3. The Feynman graphs of a few interaction processes involving zero, one, two and three condensate atoms 



FIG. 4. The Feynamn graphs of the non interaction processes 

The Green's functions can be expanded to perturbation series in the usual way. Rearranging these series one arrives 
at the generalized Dyson-Beliaev equations pO| , ^ 41 1: 



)E^/(k,zc^„)e^/(k,zu;„), 



(28) 



where is the self-energy, the contribution of those graphs which are one-particle irreducible (cannot be split to 
two by cutting a single one-particle line) and connect to two external lines with indices (r, 7) and (s,^. This equation 
is graphically represented in Fig. ^. The symmetry properties of the Green's functions of Eq. ( p^ ) stand for the 
self-energies as well. 
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y-\,\ = (-y(o)-i,-iJ(-^-i,iJ [-yii ) + (-t/(o)-i -iJ(-^-i -i)(-t/_i,ij 
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FIG. 5. The graphical symbolization of the Dyson-Beliaev equations, the hatched squares represent the one-particle irre- 
ducible graphs 



Similar rearrangements can be carried out for the perturbation series of the generalized density correlation functions 



(141). One can introduce their proper parts (the polarization parts), the contribution of those graphs (the polarization 
graphs) which cannot be split by two with cutting a single interaction line, that is they are interaction line irreducible. 
Then the equations determining these correlation functions read as: 



i?;r,, (k, ZC.„) = ffl- , (k, ZC.„) + n-,(k, lUrM'Df?,, (k, ZC.„). 



(29) 



The Feynman graph corresponding to this equation is shown in Fig. ^, where the generalized density correlation 
functions are represented as boxes which can connect to two interaction lines and their proper parts are represented 
as grey polygons also able to connect to two interaction lines. The proper parts satisfy the same symmetry properties 
as the generalized correlation functions (see Eq. (^2|) 
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FIG. 6. The graphical representation of Eq. (poj). The gray polygon represents the proper graphs 



Above the critical temperature (or for a noncondensed system), where (a,.(k)) = for all k, the anomalous 
correlation functions ( [l6| ) are zero. In this case the self-energies, which are irreducible by definition are proper as well 
and the polarization parts, which are proper by definition are irreducible as well. In the Bose condensed phase the 
appearance of the anomalous averages (ar(0)) ^ leads to the nonzero value of the anomalous correlation functions 
and to a situation where the self-energies are no longer proper and the polarization parts are no longer irreducible. 
However, they can be separated such as 



-•-yS 



^75 



M, 



■yS 7 



_ T-r(-r)sr T-r(s)sr 
— ri / / -r ri„/ „/ , 



(30a) 
(30b) 



where E is the contribution of those self-energy graphs which are proper, while the graphs contributing to M will be 
improper, and similarly II'^'"-' is the contribution of those polarization graphs which are irreducible as well and 11^'*^ is 
the contribution of the reducible polarization graphs. The separations (30a) and (30b) are the starting steps toward 
the dielectric formalism. See for the scalar gas Ref. |3l|j30|]. 



It can be directly seen from the perturbation series of the anomalous correlation functions (16a) can be decomposed 
such a way that 



A- (k,»c.„) = A-(k,zc^„)t;-(k,*c.„), 



(31) 



where is the one-particle Green's function and A^^ is the anomalous vertex, which is the sum of the irreducible 
contributions of those graphs with one incoming interaction and one incoming particle line. These anomalous vertex 
functions can be expressed with the use their proper parts and the irreducible and proper parts of the density 
correlation functions as seen in Fig. ^ 



A- (k, ZCOn) = A- (k, ZCOn) + Uif'-iK iW„)K/ A{^(k, iC.„) 



(32) 




FIG. 7. The graphical symbolization of building up of the anomalous vertex from proper parts 



C. The determination of the chemical potential and the treatment of the non-interaction self-energy vertex 

As pointed out earlier, with the introduction of the new set of creation and destruction operators with Eqs. ( ^a] ) and 
( pb| ) one can derive the relation between the condensate density and the chemical potential from the requirement that 
(6r(k)) — (6j(k)) = 0. Using perturbation theory for the evaluation of these expressions one can arrive at equations 
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which can be symboHzed as seen in Fig. ||. One can notice that a similar rearrangement is possible as was made with 
the Green's functions which leads to: 



(67(0,0)) = sg,(o,o)c;|;(o,o), 



(33) 



where Eq^ is the sum of those irreducible graphs which has only one incoming (outgoing) lines. Using the symmetry 
properties of the Green's functions one can easily derive that the consistency condition is equivalent to 



sL(o,o) = o. 



(34) 
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+ ... 



(0)l,l7 ~^^*0 ^a^a'^as ^s'V~!='(0)l 

FIG. 8. The beginning of the perturbation series of the consistency condition 



The non interaction self-energy vertex corresponding to the last term in the Hamiltonian of Eq. (|^) is to be 
calculated in each order of perturbation theory. Since this vertex is proper and diagonal in spin and anomalous 
indices it is part of the proper self-energy as illustrated in Fig. ^. 
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FIG. 9. The self-energy with the non interaction vertex 



D. Rotational symmetry and spin transfer decomposition 



Since rotational symmetry along the z axis is not broken the z component of the total angular momentum is 
conserved which results for and S!^^ in — 6s = {•y — S)(^{Fz)rsCs or symboUcally: 



, 11 ■ 



/^ii : 



py-11 • 



/->rs 

fy-ii ■ 



FIG. 10. The spin conservation property of the Green's functions. Here pQ denotes the Green's functions of the polar case 
while fQ denotes the Green's functions of the ferromagnetic case. A black box refers to a non zero element and we used 
C = (0, 1,0)^ for the polar case and ^ = (1, 0, 0)^ for the ferromagnetic case specifying the condensate 
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As a consequence of these symmetry relations Eq. (|2^) can be separated to three equations according to the number 
n = [7(0 — r) + S{c — s)]/2, with c = Cl{Pz)rsCs specifying the spin carried by the propagator relative to the spin of 
the atom in the condensate. For this reason let us define the following matrices for the polar case 



Or 

y-rS = 
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qOQ 1 




Zj^S — 
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(35c) 



where the indices in the upper left corner correspond to the number n. For the ferromagnetic one let us define the 
following quantities 



Or 

+g 


r C++ 

~ C++ 

= yi.ii 


C++ 1 
C++ 
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r V++ V++ 1 
^1,1 ^1,-1 

V++ V++ ' 
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(36a) 
(36b) 




rOQ 
= y-1-1^ 






vOO 


(36c) 


Qg 


= SlA^ 




Qe = 




(36d) 


-Qg 










(36e) 



The indices in the upper left corner correspond to the number n except for the quadrupolar spin transfer n = 2 and 
similarly —Q is for n = —2. (Note, that there is no quadrupolar Green's function for the polar case because the value 
of n can not be ±2 here.) The later modes of the ferromagnetic case (+, — , Q, —Q) differ from all the other modes 
since for these cases spin conservation forbids the existence of anomalous Green's functions and self energies. With 
these definitions one can explicitly construct the Green's functions. The propagators corresponding to different spin 
transfers will be called different modes in the following. For all three modes of the polar case {n = 0, +, — ) and for 
the first mode of the ferromagnetic case (n = 0) the Green's functions are 



— 7, — Q 



{iujn-fi ^ek - Si,i)(it^n + ^ek + S_i,_i) + E-i^iSi,-! "A 



(37) 



where ek = ek — fJ-o is introduced. This equation defines the quantities ^Na-y and "A for the first four modes. For 
the last four modes of the ferromagnetic case one can get 

±lUOn — fl Ck ~ 



±Qc - 



(38b) 



Sometimes it will be practical to cast the formulas into a most concise form. This can be achieved by introducing 
the following formal definitions: 



+ ^7^ 







7(5 








7(5 



^7(5 



^7(5 









7(5 







^-1,-1 



7(5 



i.e., the validity of Eq. 



(39a) 
(39b) 

% is extended for 



with the help of which equation (38a) takes the form of Eq. (| 
n = 0, +, — for the ferromagnetic phase as well. 

The conservation of the z component of the spin for these i^^r,/ correlation functions means that r — s = r' — s' 
which holds for the proper parts as well. And as a consequence Eq. ( p9|) also decouples into three ordinary matrix 
equations according to the specific spin transfer. For this reason it is useful to define the following matrices: 
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• For spin transfer: 



^bb 



naa 
bb 



V-aa 
bb 



n++ n++ n++ 
^++ ^00 ^ — 



n++ 

ttOO 



00 



n 



^00 

00 



n: 



i C-,-1 - Cg 

Cn Cn Cfi 
Cfi Cs Cfi Cfi -\- Cs 



where no automatic summation is understood now. From Eq. ( |29| ) one obtains 

°D = °n + "n °c °-D. 



For +1 spin transfer: 



^a)ab 


•~ ^b+l.b — 


. ^+0 


^U)ab 


-rja,a+l 
^ ^h^l.b ^ 


- ttO+ 


^a)ab 


~ ^b-^\,b — 


Cs Cs 



73)0+ 1 



^^0 



and the resulting matrix equation is: 



• For —1 spin transfer: 



B)ab 
"£)ab 



n 



a+l,a 

a+l,a 
b,6+l 



■ n+O r)+0 

^0+ ^-0 

TT + TT+O n 
yrO- ttO- 



Cs Cs 
Cs Cs 



and similarly 



• For the ±2 spin transfer case from the fact that V^^ - 



- --g -g -g, 

- V^^ = the following results can be obtained 

Knit 



(40a) 
(40b) 
(40c) 

(41) 

(42a) 
(42b) 
(42c) 

(43) 

(44a) 
(44b) 
(44c) 

(45) 



(46) 
(47) 



As a consequence of Eq. ( [2^ ) ^D_, (and their proper parts) are symmetric matrices, furthermore D_{k,iujn) — 
D_{k,~iuin) and '^D^{k,iujn) = ^Hikj—it^n)- This means that the generalized density correlation functions for zero 
spin transfer are completely symmetric under time reversal. 

Equation (p2|) also decouples according to the amount of spin transferred. 



°A 



■a 



Oj-j(^) 



(48a) 
(48b) 
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where the introduced vectors are different for the polar and for the ferromagnetic cases since the allowed spin 
projection of the incoming (outgoing) one-particle is determined by the rule of spin conservation and the spin projection 
of the condensate. It results that for a given spin transfer the spin projection of the incoming (outgoing) particle can 
take only one value (others are forbidden by spin conservation). So for the polar case (where c = Cr{^z)rsCs = 0) the 
anomalous vertex vectors are 



"A 

P—a 



A++ 

AOO 
^'-Oq 

A0+ ■ 
A+o 



A' 



0- 



.1 



-.1 



"A 

p— 



AOa 
^*-++ 
AOa 
^^00 

A""_ 



:A 



P A_, = 



aO + 



?A- 



At" 
A"- 



:a-^ = 



^^ + 

AoT' 
Air' 



(49a) 

(49b) 
(49c) 



For the ferromagnetic case (where c = C,r{-^z)rsCs = 1) it is easy to verify that spin conservation forbids any spin 
index for ^A -^ and for ^A-^. These later vectors can be taken as zero. 



fA 



+A 
/—I 



7A 
f—i 



AU\ 

A?. , 

a;«; 


0A" = 


AoY 
V A+" 


Aol ) ' 


;a^- 






7A^- 


IP 



+A 

/—-I 



tA 
/—-I 



0, 



A+O 

^*-o,-i 

-ii-n _i 



0+ 
0,-1 



(50a) 

(50b) 
(50c) 



for the ferromagnetic case. The symmetry relations of Eq. ( ps] ) hold for the proper and irreducible parts as well, 
which is equivalent to: 



°A"(k,zc.„) 
+A"(k,iw„) = 



°A"(fc,ic^„) = "AJfc,iL^„) = °A_Jk,-iLu„) = OA "(fc,-jc^„), 
+A"(fc,ici;„) = +A^(fc,«w„) = "A_^(fc,-ici;„) = "A "(fc,-iw„). 



(51a) 
(51b) 



E. Dielectric functions 



With the definition of the 



dielectric functions Eqs. (^91), (41) 



°i(k, zc^„) = "i-°n(k,ic^„) -Ofi, 

^e(k, iw„) = ±1 - ±n(k, iujn) ■ 



}43| ) and ( [4 51 ) can be rewritten as 

e^;;(k, iuJn)D% (k, itj„) = fi n^r,, (k, iu;„), 
°£(k, iw„) • °0k, iw„) = On(k, zw„), 
=^e(k, iw„) • =^i2(k, it^„) = '^n(k, zcj„). 



(52a) 
(52b) 
(52c) 

(53a) 
(53b) 
(53c) 



F. Interaction propagator 

With the use of the proper parts of the density correlation functions one can define an interaction propagator: 



' y r' s' * r' s' ' ' ^ ao cd* r' s' 



(54) 
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which can be symboHzed as depicted in Fig. ^ Using the dielectric functions one can get 

>V:,^eX, =Kv, (55a) 
This equation also splits to parts according to the spin transfer. 

0^ 0| = 0£, (55b) 

±^ ±£ = (55c) 

Wl+ = = 0. (55d) 



+ <: 



r s 



r s 



r s 



FIG. 11. The Feynman graph of the interaction propagator 



G. Proper Green's function and irreducible polarization part 



With a simple substitution one can verify that both Eg. (|28| ) and Eq. (55a) have the following partition property. 
If one decomposes the self-energy and the proper part of ( p_41D into two 

^^s-^T + ^T' (56) 

n^r,, = + ^(?lf^ (57) 

then one can define propagators corresponding e.g. to the first part of these quantities 

w^l:- = v;,i, + wTuTVrt' (59) 

in such a way that 

, = will:' + w:tu^;'j'"'wi}if' (ei) 



is fulfilled at the same time. Specially if one decomposes the self-energies by ( pOaD and the proper graphs by ( pOq ) it 
defines the proper Green's functions 

and an effective potential 

w,7^, = v;,i, + w^^uif'^vX- (63) 

such a way that 



g;t^g;i + g;iAC;g;f, (64a) 



= W,7,, + W:tIlif''W,% (64b) 
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is fulfilled. The earlier discussed symmetry properties hold for these decom pose d parts as well (since these are defined 
by a class of graphs) and as consequence, Eqs. (|6^), ( |64a|) and Eqs. (|63|), ( |64b| ) split to matrix equations in the same 
way as Eq. (g8|) and Eq. (|2^) did. It is convenient to define the regular part of the dielectric functions with the 

- m - ^TVc'a (65a) 
o/)^o^_OgW.o^^ (65b) 

±/)^±l_±nW.±^ (65c) 



equations which can be used to express the effective potential as: 



(66a) 

(66b) 
(66c) 



H. Improper self-energy and singular polarization 



With the use of the irreducible and proper anomalous vertex functions and the effective potential one can construct 
the improper self energies (as seen in Fig. p^.): 

n M^^(k,ic^„) = A''^20^,iu;^)W^f{-k,iujr.)Ai;{k,tujn), (67a) 

ft "M^p = °r 0^ % = or °A;3, (67b) 

ft ±M„^ = ±r ±^ ±A^ = ±r ±£ ±A^. (67c) 

It is easy to verify that = which means that for the ferromagnetic case the improper self-energy in the ±2 
(quadrupolar) spin transfer mode 

^^M = 0. (67d) 



m 





b a 

/3 " 




c 9k 



d hi f 




+ 



FIG. 12. The structure of the perturbation series of the improper self-energy 



The singular polarization can be expressed as well with the use of the proper Green's functions and the proper 
anomalous vertices as seen in Fig. O. and is of the form: 



ftn^f^f (k,»c.„) = A- (k,zc.„)CMk,*c.„)A^f,,(k,zc^„), 
n^R^'^^ ±^„,±A„o±A^ 



(68a) 
(68b) 

(68c) 



where the circle denotes the diadic product operation. For the polar case the singular polarization in the ±2 spin 
transfer mode is equal to zero since there is no corresponding proper Green's function. For the ferromagnetic case 
e.g. for the +2 spin transfer mode the singular polarization is: 



ft fU 



(69) 
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FIG. 13. The diagrammatic structure of the singular polarization 



As for the scalar case one can see intermediate states corresponding to collective excitations in the perturba- 

tion series of the one-particle propagators and at the same time one can identify one-particle intermediate states in the 
perturbation series of the correlation functions of the collective modes. However there is a great difference between the 
scalar and the spinor models. The rotational invariance of the interaction potential results in that = = 
which leads to the proper nature of the quadrupolar spin density correlation functions. This combined with the fact 
that in the polar case there is no anomalous vertex with ±2 spin transfer results in the prohibition of the coupling to 
any of the one-particle correlation functions in this phase. 



I. Couplings amongst the correlation functions 



With the help of the quantities discussed in the previous subsection the properties of the density correlation 
fimctions can be further investigated. For this reason let us introduce the following vectors: = (1,1,1)VV3, 
= (l,0,-l)V%/2, ^3 = (1,-2, 1)VV6 for the zero spin transfer mode and = (1, 1)V^ and = (1,-1)VV2 
for the ±1 spin transfer mode. It is easy to verify that {Cj^iCj'^J" ^"^^ KX-^tX^ are orthogonal and normalized 
basis s ets i n the linear vec tor spaces with dimension three and two, respectively. The density correlation functions 
(14a),([l4b|),(14c) and (14d) can be expressed as follows: 



C„„ - 3 (el , i^.. - 2 °£ , = V6 (^1 °^ , 



(70a) 
(70b) 



Since "Z? is a symmetric matrix Dnz = D^n holds. 

Multiplying Eq. ( ^ ) with 2%^ from both sides and using that = 2csX^ ° Xi oi^s obtains 



fin. 



'-n++D++, 



(71) 



where 11++ = 4xJ^nXj^. A similar equation can be derived to D There are both regular and singular contributions 

to these propagators. 



One can get a closed system of equations for the zero spin transfer correlation functions of Eq. (14a), Eq. ( |l4b| ) 
and Eq. (14c). To this end first we multiply Eq. (41) with \/3Cj^ from both sides. We further use that — 

(72a) 



Scn^j^ ° ^1 ~'~ ^"^"^2 ° ^2' "^^^ resulting equations sound as 



where n„„ = 3^^°!!^^ and Ilnz — V^S.t^IH^- Multiplying ( p] ) with v^^^ from both sides one gets 

Dzz = flUzz + CnllnzDnz + CsUzzDzz, 



(72b) 



with 1X^2 = '^^^'^JM,^- The third equation can be obtained by multiplying Eq. ( pl| ) with y/S^^ from the left hand side 



and with \/2(_^ from the right hand side which leads to 



(72c) 



These coupled equations can be solved for the correlation functions. The polarization part Ilnz written out in detail 
reads as 



n 



n 



It 



00 



n' 



00 



(73) 
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It is obviously zero if the system is invariant under spin reflection, i.e., when those matrix elements coincide which 
can be get from each other by reverting a + to a — and vice versa. This condition is fulfilled for both cases in the 
s ymm etric pha se, moreover for the polar case it even holds throughout the condensed phase. With Unz — Eqs. 
(72a) and ( 72b ) are independent, giving the solutions 



Dnr. 



fin. 



1 

0, 



fin. 



1 - Csli^z 

and leading to two separate excitation spectra. The singular part of Yizz can be cast with Eq. (|E 

fin(«j = 2 (e^"A J («A^e,) 



(74a) 
(74b) 
(74c) 



]) to the form 



(75) 



which is zero if spin reflection symmetry is present, since the anomalous vertex vector ^A^ is orthogonal to this 
case. 

In the ferromagnetic phase the spectra of collective excitations corresponding to spin density waves and density waves 
with n = are coupled. Above the critical temperature this coupling vanishes and the two modes will be independent. 
In the polar phase spin density fluctuations will always be independent of the particle density fluctuations, furthermore 
lii'} = which means that the spin density correlation function (|74g) has only regular contribution. 

In the polar phase further analysis leads to results important for the following. First note that, when the system 
exhibits the spin reflection symmetry, '^Z? s^nd ''n have only 4 independent elements (instead of 6 characteristic to a 
symmetric matrix). These independent elements are ntt = nil. nVT = nil, n1°, = li^^ = nH = nQQ^ and 

ttOO 



Second, these can be cast into a more practical form with the help of the orthogonal matrix 



h 

2 
6 



v2 

3 



3 



k 

2 

6 



(76) 



Its rows are built from the vectors {^j^: ^s-"^' Performing the transformation leads to 



Oj-[(a:)' 



3c„ 
2cs 




J0 . u^v 














3 


3\/2 







2 












^^33 




- 3\/2 


6 J 



(77a) 



(77b) 



t(^) 



with nia^ = 3V2C^°n ig and li^^ 

parts, respectively. Furthermore, according to 
The total polarization matrix therefore becomes 



6^g'^n*'^''^g, where x can be r or s corresponding to the regular or the singular 



I) ni^i - 0. 



V = 



3 



3\/2 



2 



3%/2 





(78) 



Using Eqs. (77a) and (78) the dielectric function (52b), after the transformation, can be cast into the form 



- c„ (n,^„^ + n 





c.ni';^ 



t('-) 
n3 



n 



(79) 
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We arrive at the important result, that the determinant of the dielectric function, which is invariant under such 
transformations, reads as 



det°e = det°£' = 



1 — c. 



(80a) 



This determinant factorize into two in agreement with the separation of the density and spin density fluctuations. 
Furthermore it shares one factor with its regular counterpart since 



det^e^'^^ 



1 



(80b) 



which can be obtained by taking the determinant of Eq. ( [79|) after setting all II'*^ = 0. 

Another i mpor tant consequence of the spin reflection symmetry i s th at the interaction propagator (55b) and its 
regular part (66a) has the same structure as the interaction matrix (77a): 



3Cn 
2Cs 











d'-^ 





(81a) 



with 



Cn 
Cs 



1 - CnUr, 



1 - r n*''^ 



1 - c.ni'J ■ 



(81b) 
(81c) 



Those correlation functions which have singular polarization parts correspond to fluctuations belonging both to the 
condensate and to the noncondensate, while those which have no singular proper part only belong to the noncondensate 
(they can not couple to any of the one-particle modes). So in the polar phase the spin density fluctuations described by 
Dzz (which are independent from the particle density fluctuations described by £*„„) belong only to the noncondensate. 
The singular polarization parts of the quadrupolar modes are also zero (since there is no corresponding Green's function 
and the anomalous vertices of the ±2 spin transfer are also zero) meaning that these type of fluctuations also belong 

only to the noncondensate. The proper parts n„„, n++ and 11 have singular contributions too with intermediate 

one-particle states ^Gai3, ^Gap a,nd ~Qaf3, respectively. These modes are coupled to the condensate. In the 
ferromagnetic phase Dzz and Dnn a,re coupled (their cross correlation Dnz do not vanish either), their denominators 
are common. All of their proper parts have singular contributions resulting in their coupling to the one-particle 
excitations with ^Q. These fluctuations thus belong both to the condensate and noncondensate. The correlation 
functions of the spin waves (both the ±1 and ±2 spin transfer modes) Z?--, -DJ_|_ and D?_ have singular 

proper parts connecting them with the propagators ^Q, ^G, '^G a-nd ^'°^G, respectively. Note, that there are no 
improper self-energies for the quadrupolar spin transfer modes, consequently their proper Green's functions are the 
same as the full propagators of these two modes. 



J. Coupling between the Green's functions and the correlation functions and the spectra of excitations 



A) Let us first consider the modes with n = 0, — . With the decomposition ( pOal ) of the self-energies and with a 
straightforward calculation starting from Eq. (5^) one arrives at the expression 



'A "A - "iV^r "Mt^ - det "M' 



Ga-r = „ ;v = ~T ~ — T7Z , , „ . , ^ (f^^a) 



both for the polar and ferromagnetic phases. Here we introduced the 

"TVa^ = 5a-y{aiujr, + fi^^ek) + a-f (82b) 
"A = (ic^„ - fi-^ek - "Sn)(ic^„ + fi^'ek + + (82c) 
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quantities, with = "7V„^/ " A and det "M = "Mn "M_i,_i- "Mi,_i "M_i,i. With the help of Eqs. ( pb| ) 

and ( 53c ) the density correlation functions with n = 0, +, — can be brought to the form 

= = „= , (83a) 



by writing 



= det "£ 



^S = det("e) 



(83b) 



An important relationship can be shown in the following way. First applying the decomposition "e = "e*-'"'' + "e*-^\ 



with = to det "i leads to 

det "£ = det ("l^'^^) det (l - "g^"^ 



(84) 



With the help of Eqs. (B8t) and (38c) and a straightforward but rather lengthy algebraic manipulation one can find 
that 



det ( 1 - "n^') "M^ ) = det 



i - fi"' "Sa/J^A^ o "A'' = 1 - fi-i "ea/3"A„ "A' 



fi-i det ( "^g) det ("A^ "A ) . (85) 



Inserting it to Eq. ( ^4) ) and using Eqs. ( |67b|) and (67c) and that det — — "A ^ and comparing the result with 
the denominator of the Green's functions ( ^2a|) one can arrive at the basic connection 



'Adet"e 



^Adet 



(86) 



Changing to retarded correlation functions, which is done in the usual way (by analytically continuing in fre- 
quency) pO[ , the elementary excitations of the system are determined by the poles of the corresponding correlation 
functions. The spectra of the one-particle excitations (quasiparticles) are given by the poles of the Green's functions, 
or equivalently by the equation 



'A = 0. 



(87) 



The spectra of collective excitations are determined by the poles of the density correlation functions, or equivalently 
by the equation 

det = 0. (88) 

Equation ( |8^ ) shows that the Green's functions with spin transfer n = 0, +, — can be arranged to have the same 
denominator as the density correlation functions with the same spin transfer. This means that if det "e has a zero (the 
corresponding density correlation function, "Z?, has a pole) than "A must have a zero ( '^Q must have a pole) there 
as well in general. The zero spin transfer mode in the polar case is an exception, where both the dielectric function 
and its regular part factorize and share some of their zeroes (see Eqs. (SOa) and (pOq)). In this case the zero of the 



left hand side coming from det satisfies the equation with the zero coming from det '^e^'^^ of the right hand side, 
instead of °A. Hence, this zero do not appear amongst the poles of ^Q. This further means that density correlation 
functions with and ±1 spin transfer has common denominators with the Green's functions corresponding to the 
same spin transfer for both the polar and the ferromagnetic cases, except the Dzz spin density correlation function 
in the polar case, which do not couple to any of the Green's functions. 

B) Concerning the quadrupolar spin waves (n = ±2 modes) let us consider first the ferromagnetic phase. Both the 
full Green's function and the correlation function are proper (see Eq. (|67d|) , (^ ) and (|4^)) for such spin transfers. 
Consequently one can subtitutc in Eq. (|69|) the full Green's function leading to 



n 



(r)T± 



,±1 



(89) 



HID therefore these type of collective 



The one-particle excitations then appear also as collective ones. 

In the polar case there is no Green's function with n — ±2, see Sec. 
excitations, do not show up as one-particle elementary excitations. 

Thus in the polar case the collective and one-particle excitation spectra do not match fully. There are excitation 
modes in this phase, which do not belong to the condensate. In the ferromagnetic case all of the collective excitation 
modes are connected to corresponding Green's functions, meaning that all types of density fluctuations are governed 
by condensate dynamics. 
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IV. BOGOLIUBOV APPROXIMATION 



In this section we discuss the simplest approximation in the framework of the dielectric formalism, namely the 
Bogoliubov approximation. This approximation is a very low temperature calculation, it neglects all terms coming from 
the non-condensate density. This approximation was studied by other authors as well with other techniques. 

With the Green's function method one can give the form of the correlation functions as well not just the frequency of 
the modes. 

The proper self-energies in the Bogoliubov approximation are 



The corresponding Feynman-graphs can be seen in Fig. ul 



(90a) 
(90b) 



a 



+ 



r' Cps' 



r I 



+ 



s' Cp r' 



sO- 



FIG. 14. The Feynman graphs of the proper self-energies and the tadpole diagrams in the Bogoliubov approximation 
The regular polarization is zero in the Bogoliubov approximation and the proper anomalous vertex reads as 

= ^ [SraSc-lCs + 6saSc.,lC] , (91) 

which can be graphically represented as seen in Fig. n5l 




+ 



s a 



FIG. 15. The proper anomalous vertex of the Bogoliubov approximation 



After calculating the prope r Gr een's functions the improper self-energy and the singular polarization can be con- 
structed from Eqs. (67a) and ( |68a| ) . With their explicit form known the correlation functions can be expressed. These 
calculations have to be made separately for the polar and for the ferromagnetic cases. 



A. Polar case 



In the polar case Cc, > 0; the condensate spinor in this phase can be chosen to = (5r,o- The chemical potential 
can be calculated using Eq. (^ which results in ^ = N^Cn. Since the proper self-energy is zero (except the diagonal 
fiQ terms) the proper Green's functions will be: 
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This is similar to the free propagator of Eq. ( |27| ) but the /io parameter has been canceled out. The anomalous vertex 
vectors are: 



OA 

P—o 




?4 











(93) 



With the use of them the singular polarization matrices and the improper self-energies can be expressed using Eqs. 
(pq) and (|6^). The singular polarization matrices are 




Hs 








n± 



(94) 



where Il± — No/{±ifkUn — £k) and Us = 11+ + n_ were introduced. The improper self-energy matrices are: 



NoCn NoC„ 
NoCn NoC„ 



= r 



NoCs NqCs 
NqCs NqCs 



Putting all together to Eq. (B2a) the Green's functions will be 

Sa-yiaiujn + h^^Ck) + ajh^^ NoCn 
" (*a;„)2 - fi-2ek(ek + 27Voc„) ' 



y. 



7 



(icj„)2 - H 2ek(ek -I- 2NoCs) 



(95) 

(96a) 
(96b) 



By the use of Eqs. (g3|) the correlation functions (40a), (42a) and ( |44a| ) can be cast to a form, after multiplying both 
the numerator and the denominator with "A — {iun — fi~^ek)(i'^n -I- fi^^ek), where the denominators are common 
with the corresponding Green's functions in agreement with Eq. (|8^)) 



fi-i27Voek 



„)2 _ fi-2ek(ek + 2iVoc„) 
1 

— (icj„)2 - /i-2ek(ek + 2iVoCs) 




1 




-D = 



(97a) 



(97b) 



With the use of Eqs. IxM) the particle number and spin density correlation functions can also be calculated. The 
result is that D„ti = 2fi A^oek/[(«w„)2 — fi~'^ek{ek + 2NQCn)] while D^z = Dm = 0. For spin waves D++ = = 
4fi-iiVoek/[(ia;„)2 - H~'^ek{ek + 2NoCs)]. From Eq. (M) and Eq. (|68a|) d2± = as mentioned before. 



B. Ferromagnetic case 

In the ferromagnetic case Cs < 0. One can choose the condensate spinor as = Sr,+- The calculation is analogous 

to that made in the previous subsection. The chemical potential (from Eq. (|^) is ^ = No{cn + Cg) = Nog. The 
proper self-energy has 3 different components in its diagonal (other self-energy components are zero) which results in 
the following proper Green's functions: 



aiujn — fi ^Ck 

C^(k,^c.„) = — (98b) 

So 

„-fi~Hek-2iVoc,)' 



aiLUn - fi^^(ek - NoCs) ' 
g~~ (k, tUn) = — ^ . (98c 



There are two things to note. The first one is that the parameter has cancelled out here as well, the second 
one is that the remaining self-energies (in the last two proper Green's functions) are negative since Cs < in the 
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ferromagnetic case. The anomalous vertex vectors are: 




+A 
f—i 



The singular polarization matrices are: 



on(«) 



nil++ 








n 



is) 



n: 







+0 





(99) 



n 



n(^+° 





with 



n 



{s)++ 



1^ 



n: 



(s)+0 
0+ 



The singular polarization for the ±2 spin transfer reads as: 



n 



The improper self-energies are: 



7Vo5 Nog 



(100a) 

(100b) 
(100c) 
(lOOd) 

(lOOe) 
(101) 



where </ = c„ + Cg. For the ferromagnetic case the resulting Green's functions are (for the later ones see Eqs. 



°G - a++ 



yii = yii = - 

I' 



1 



1 



iLUn - h i(ek - 2NoCs) ' 



(102a) 
(102b) 
(102c) 



For spin transfers n = 0,+,— the correlatio n fun ctions (40a), (42a) and (44a) can be evaluated from Eqs. 
and with the use of the polarization matrices (100). Both the numerators and denominators are multiplied with the 
corresponding "A, leading to 



3^ ^ 2fi-^iVoek 

= ~ (zw„)2 ~ n-^euiei, + 2Nag) 

No 



1 





D 



iojn — fi 
No 

—iuiji — fi~^ek 



1 




1 




(103a) 

(103b) 

(103c) 
(103d) 



The correlation functions with n = ±2 obtained directly from Eqs. (^6|) and (|4^) with polarization functions (lOOe 
read as 



(103e) 
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The particle number and spin density correlation functions from Eqs. (uG) are: 



[iLUnY - n ^ek(ek + 2Nog) 
^++ = -^^, (104b) 



D-= . , (104c) 

-lUJn - n ^Ck 



C. Collective excitations in the Bogoliubov approximation 



The spectra of collective excitations can be expressed using Eq. (jS^) (for the and ±1 spin transfer modes) or 
equivalently by the zeroes of the denominators of the appropriate retarded correlation functions. For the polar case 
this results in 



"iu = fi-i Vek(ek + 2Noc„) ^ J (105a) 

V 



M 

^iu = fi-Vek(ek + 27VoC.)^-^^ V^^^' ^^^^^^ 

There are no other excitation modes for the polar case in this approximation. These modes have linear spectrum and 
are Goldstone modes. The °w mode belongs to the '^Gaj Green's functions and the -D„„ (^D.) correlation functions, 
while the ^ur mode belongs to the ^Ga-y Green's functions and the D±± (^-D) correlation functions. 
For the ferromagnetic case the excitational energies are at: 



°cj = fi" Vek(ek + 2iVo5) ^ \l (106a) 

= ift-^ck, (106b) 
^Lj = ±n~^{ei,-2Noc,). (106c) 

The first {'^uj) mode, which is responsible to density (and spin density) fluctuations, belongs to the '^Ga-i Green's 
functions and the connected Dnn, D^z and Dnz CM) correlation functions. This is a linear Goldstone mode. The 
next {^iv) modes are responsible for ordinary spin waves and belong to the ^Gii = G±^± Green's functions and the 
D±± (^D) correlation functions. These are also Goldstone modes but with a quadratic dispersion. The last (gw) 

modes describe quadrupolar spin waves. They belong to the qGii = G±± Green's functions and the -D^i (^±t) 
correlation functions. These are non Goldstone modes and they start with a gap. The frequencies found in the 
Bogoliubov approximation agree with those of Refs. lUJl^ obtained in a different manner. 



V. RANDOM PHASE APPROXIMATION 



The simplest way to take into account the appearance of the noncondensed atoms, leading to a damping mechanism, 
is to add a Hartree term to the Bogoliubov proper self-energies (compare Fig. ^ with Fig. |l^) and to choose the 
regular polarization graphs as bubbles. Some results of this approximation were published earlier in pl[ |. 

Therefore the proper self-energies of the random phase approximation can be graphically represented as seen in 
Fig. 



16 and their contributions are 



(107a) 
(107b) 
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Here the i?'"' notation is introduced for the contribution of the Hartree term: 



r s 



r s 

-Ar- 



+ 



r'qps' 



r I s 
I 



r' Cp s' 



+ 



r I s 
I 



+ 



IS 



+ 



IS 



r O- 



r O- 



FIG. 16. The graphical representation of the proper self-energies and the tadpole graphs in the RPA 



The proper self-energies are independent of the wave-numbers and frequencies and diagonal both in the Roman and 
Greek indices. The proper Green's functions are diagonal in the Roman and Greek indices as well. The approximation 
is done in a self-consistent way, namely the internal propagators used are proper Green's functions in the Hartree 
approximation. 

The regular polarization will be a bubble (with proper Green's functions) as seen in Fig. |l^. Its contribution is: 
-nU^;,lr{k, zu;„) = J E '(q- *^")^n (k + q, ic.„ + *i.„). (109) 



Since our approximate proper Green's functions are diagonal in their spin indices (r = r' and s = s') the regular 
polarization matrices will be diagonal as well, with 



H 



('■)++ 

f+ 






h: 





(r)00 
00 





H 








H 



(r) 



H 



'r)0+ 
+0 





H, 





(r)-0 



(110) 




FIG. 17. The Feynamn graph of the regular polarization function in the random phase (Hartree) approximation 

The anomalous vertex remains the same as it was in the Bogoliubov approximation, see Fig ^ and Eq. (pl|). With 
the building blocks determined we can now turn our attention to the construction of the correlation functions. 



A. Polar case 

Using the interaction (Q) with Cs > and Cr = Sr.o the self-energies (107) in the polar case 

E++ = h-' [fio + {H++ - H—) Cs] , 
^11 — A*0, 

Err = h - {H++ - H—) . 



(Ilia) 
(111b) 
(111c) 



Since in the polar case the system has zero magnetisation the Hartree terms should satisfy = H and as a 

consequence the proper Green's functions are the same as they were in the Bogoliubov approximation Eq. (|9^). The 
situation is similar to that in the scalar case, the chemical potential cancels the Hartree terms |pl|. 
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The anomalous vertex vectors are given by Eq. (^), and so the singular polarization matrices will also be the same 
as for the Bogoliubov approximation Eq. (|94|). The regular polarization matrices are as follows: 



Ho 
Ho 
Ho 



n 



(r) 



Ho 
Ho 



in accordance with the spin reflection symmetry discussed in Sec. [III. Here the notation Hq = 11' 
'"''^ is introduced. The interaction propagators are then obtained from Eqs. (66a) and (66b), 



(■'■)++ 



n 



(112) 



(r)00 



n 



(r) 



ri^) 



(r) 



,(r) 



_ nV'l n\i') nV 



{r) p{r) p(r) 



w = 



with 



CW(k,zc.„)- 



1 - 3c„no(k, ?w„ 
With Eq. (|67|) the improper self-energies become 



CM(k,. 



'-M = h- 



1 - 2CsUo{'k,iu!n) 



Nod'-'' Nod"^^ 

Nod"-^ Nod"-^ 



(113a) 



(113b) 



(114) 



These can be used in Eqs. ( |82| ) to arrive at the Green's functions in RPA. After multiplying both the numerator and 
denominator with det "e^''^ we obtain 



Or 

ya-) 



6aj{aiLUn + fi ^ek)(l - 3e„no) + a7fi ^NpCn 
[{lUJn)^ - fi-2e2](l - 3c„no) - 2fi-2ekiVoc„ ' 
6aj{aiL0n + H~^ek){l - 2csno) + ajh^^NoCg 



(115a) 
(115b) 



Calculating the regular dielectric functions from Eqs. (|65|) and using Eqs. (|83| ) the correlation function matrices 
could be explicitly given. But in the polar case one can use Eqs. ( [t^ to calculate directly the relevant correlation 
functions in the zero spin transfer, with n„„ — 3c„no + Cnlls and IIz^ = 20^110 leading to 



[{iujnY - fi-^el]{l - 3c„no) - 2fi-2ekiVoc„ ' 

2no 
i-2c,no' 



Dnz = 0. 

With Eq. ( [7l|) and with n±± — AWq + 2ns the correlation function of the spin waves will become 

4no[(ic^„)2 - h-^el] + 2H-^Noek 



D±± = h— 



-2^2 



](l-2csno)-2fi-2ek^oc/ 



The correlation functions of the quadrupolar spin waves come from Eq. (p9f) with a result of 



D 



Q 

±± 



MUq. 



(116a) 

(116b) 
(116c) 



(116d) 



(116e) 



One can easily verify that the denominators of the appropriate correlation functions match as was shown at the 
general formalism. 
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B. Ferromagnetic case 



In this phase Cs < and Cr = Sr,+ and with the interaction the self-energies (107) become 

= fi-i (/io - Ale,) , 



where the newly introduced quantity 



M=No + H++ - H- 



(117a) 
(117b) 
(117c) 

(117d) 



is the total magnetisation of the system. It describes the number of particles responsible for the magnetic mean-field. 
Note that since the Hartree terms are momentum independent the proper self-energies will also be. Therefore the 
proper Green's functions for the ferromagnetic case read as 



aiu}„ 



fi-i(ek - 2Mcs)' 



(118a) 
(118b) 
(118c) 



These proper Green's functions are similar to their Bogoliubov counterparts only the momentum independent mean 
field has changed with the appearance of the noncondensed particles with spin projection -I- and — . The chemical 
potential here do not fully cancels the proper self-energies leaving behind the energy shift due to the magnetic mean 
field. 

The anomalous vertex vectors are given by Eqs. (|99|). One can express the singular polarization matrices from Eqs. 
( 3^) f or n = 0, -f,— , and they are of the same form as their Bogoliubov counterparts (100a) with elements ( |100b[ ), 
( 100c ) and (lOOd) where the proper Gree n's fun ctions (US) are to be used. The singular polarization for the ±2 spin 
transfer can also be brought to the form (lOOe) with (118c) as the proper Green's function, 



(119) 



Using Eqs. (|65|) and Eq. (110) the effective potential W can be constructed. Inserting it to Eqs. (|67|) the improper 
self-energies will read as 



= detOe^"' 



1 1 
1 1 



1 




dct 

where the following notations are introduced to simplify the equations 



0, 



(r) — 



n, 



(r)OOs 



(i - c„nM) (i 



c c nw 



dei+e'" = 1- 



CS yr(r) 

2 ++' 



with 



n(-) 

nn 


= n 


r)++ , 
f+ + 


^(r)00 „(r) 


n('-) 


= n 


r)++ , 
f+ + 




n(-) 

nz 


= n 


r)++ _ 




n(-) 


= 2 


[kr 





(120a) 



(120b) 
(120c) 
(120d) 

(120e) 
(120f) 
(120g) 
(120h) 
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From Eqs. 



and from Eq. (38b) the Green's functions are 



Or 

+^11 



- r+^ 



^00 

= yii = 



[(iw„)2 - fi-2e2 ] det °e^''^ - 2n-^euNop ' 



det+e 



^Gii = Gil 



fi-i (ek - Mc,)] det - fi-iiVoCs 
1 



(121a) 

(121b) 
(121c) 



Equations (33a) and (|89| ) give the correlation function matrices which combined with Eqs. ([70|) result in 



Dnn = ^ 



n 



++ 



1 - c, n 



.(r)00 



4n 



(r) — 



n, 



(r)00 



n 



(r) — 



1 - c^Hoo 



(r)00 



'n 



ir)++ 
++ 



n 



+^ 



det 

i-c„(n(r + 4nL^^" 



) [l - c„ (l 



n' 



(r) — 



1 - c n^''^"° 



det^e 



n 



('■)++ 
++ 



det^e 



D'l 



n 



(r-)0+ 



2fi- 



n 



(s)0+ 
+0 



(r-)+0 
0+ 



det+£ 



(122a) 

(122b) 

(122c) 

(122d) 
(122e) 

(122f) 
(122g) 

Equations 

(122a), (122b), (122c) and (122d) can be cast to a form having common denominators with the corresponding Green's 
functions by multiplying both their numerators and denominators with the apropriate "A. Only writig down the 
resulting denominators: 



Ah It 



Tir)- + 



n 



+- 



Here 



det °£ = (1 - C„n„„) (1 - CsU^z) ~ CnCsli. 



2 

nz J 



det +e = 1 - 7fn++, 



where n„», = 11' 



(r) 



n 



s)++ 

++ 



, n,, = n 



(r) 



tt(«)++ n 



nil + n^|++ and n++ = nit 



n: 



(s)o^ 



° A det °e = [{iujn f - h-^el] det - 2h-^ev,Nop 



"A det +e = [(iw„) - fi"^ (ck - Mcs)] det +e 



(r) 



(123a) 
(123b) 



For the rest of the paper we will deal with the retarded correlation functions which can be obtained in the usual 
way by analytically continuing in frequencies (see e.g. Ref. [Q). First however we will discuss shortly the static 
properties. 



C. Static properties of the spin-1 Bose gas in the random phase approximation 



We will study the equation of state and connect it to the density autocorrelation function through the compressibility 
sum rule. The equation of state is investigated by choosing p, N as conjugated variables. The number of particles N 
is given by 

N^Na + H"", (124) 
where iJ'*'" is defined by Eq. (108). In the ordered phase Eqs. (^), ( 117(^ ) and (124) provide a second relationship 

p = CnN + CsM. (125) 
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It is convenient to introduce the following quantities: the thermal wavelength A = fi(/3/2M)-'^/^, the critical particle 
density of the ideal gas Nc = ST (|) C(i)/[(27r)^A'^], dimensionless interaction strengths en,s = NcP\cn,s\, relative total 
magnetisation m = A4 /Nc, relative particle number x = N/Nc and dimensionless chemical potential u = (3fj,. Here 
r(s) is the Gamma function and C(s) is the Riemann-zeta function. 

Both for the polar and ferromagnetic cases in the symmetric phase (where no condensate is present and the total 
magnetisation is zero) the equation of state reads as 



1 



C(i) 



i^(|,a;e„-u) 



where 



1 



oo J-S — 1 



dt 



is the Bosc-Einstein integral. 

In the condensed phase for the polar case the chemical potential is determined by Eq. 



|) which leads to 



CnX. 



(126) 



(127) 



(128) 



Equations (126) and (128) can be used to give the isotherms u{x) of the polar case. Equation (126) is valid when 
X < 1 and Eq. (128) is to be used if x > 1. These isotherms show the character of a continuous phase transition. 
In the condensed phase for the ferromagnetic case the chemical potential can be calculated also from Eq. (p3) 



sm{x) 



i{x)], 



(129a) 



where k = Sn/^s and m{x) is the relative total magnetisation as a function of the relative total density given by the 



3C(i) 



,) + 2F(|,2me,)] 



(129b) 



equation. The isotherms are given by Eqs. ( 126 ), (129a) and ( 129b|) . The solution of Eq. ( 129b ) is unique if a; > 1 
but at a; < 1 a second solution emerges also with nonzero magnetization and both solutions can be continued down to 
Xy < 1, where they coincide and vanish. This means that the chemical potential has three solutions between x — Xy 
and X = 1. (One of them with zero magnetization at a given density.) In the applications Cs ^ 1, so the Bose-Einstein 
integral can be well approximated by i^(|, a) « C(l) ~ l^pKot leading to 



At X = 1 the two solutions are at 



= 1 - 



'^1,2(2^ = 1) 



(1 + 2\/2)27res 



9C( 



3\2 



4(l+2V2)^7res 

0. 



(130) 



(131) 



These isotherms show the character of a very weak first order phase transition. A similar equation of state with a first 
order transition was derived in p9| for scalar condensates but in a different model approximation. In that work the 
first order character of the phase transition was caused by the exchange interaction (Fock term) which is known to 
lead to such behaviour The situation is interesting, since for scalar particles the RPA (Hartree approximation) 
gave a continuous phase transition as here for the polar case 30|. 

Since the above behaviour of the equation of state is not common a further consistency check is required. The 
compressibility sum rule is a candidate, since it relates the static correlation functions to the derivatives of the 
equation of state ||4^,|33| , namely 



=ft-iAm(fc^O,c^ = 0). 

0[i J T 



The inverse of the left hand side of Eq. (132) can be easily calculated, since 

m'(a;)es] . 



dpi 
'dN 



1 /, N 1 

u [xj — -z-rr iKC.s 



(132) 



(133) 
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The derivative m'{x) can be expressed from Eq. (129b) with the identity daF{s, a) = —F{s — 1, a), leading to 





"<o)"°(o,o) + 4n(:l-(o,o)" 


1- Cs 




^°"(0,0) + 40^1"" (0,0)' 





(134) 



where we used that n[,o^''°(0, 0) = -/3A^c-F(i, TOe^)/(3C(i)) and U^J^l (0,0) = -PNcF{^,2mes)/(3C(^)). It is 
straightforward to verify that the static hmit of H/Dnn is just the right hand side of Eq. (134) in agreement with Eq. 
( |132D . The point Xy where the two solutions mi and m2 coincide is the point where the graph of u{x) has a vertical 
tangent. The point where the graph of u(x) has a horizontal tangent is at where the density autocorrelation function 
has its pole at g = and a; = 0. This happens at 



Xh = l- 



2k (1 + 2^/2)27res 



(«:-l)2 9C(i)2 ■ 

For X > Xfi (in the condensed phase) the static correlation functions have no real positive pole in the q variable. 



(135) 



D. Collective excitations in the random phase approximation 

The excitation spectra can be obtained in the usual way as the poles of the retarded Green's functions and correlation 
functions. The real parts of the poles correspond to the energies while the imaginary parts correspond to the damping 
(the inverse of their lifetimes) of the excitations. 

Our further calculations can be done more appropriately in a dimensionless form. For this reason let us further 
introduce the following characteristic lengths: ^ = fi{4:AI No\cn.s\)^^^^ the mean field correlation lengths associated 
with the interaction strength c„ and Cg, and ^' = M / {AttH^Nq(3) the length scale of the critical fluctuations. The 
calculation will be limited to the intermediate temperature region. This region is defined by the ^ 3> A conditions. 
In this temperature range the fcA <C 1 condition can be fulfilled for the physically interesting wave numbers and the 
contribution of the bubble diagrams can be taken as perturbation for the frequencies considered. 

A dimensionless frequency can be introduced with 



= ^. 
ek 



The singular polarization functions in the polar case take the following form: 

c„,sns(k, w) 



CB 



1 



1 



(fcA)2 



1 



While in the ferromagnetic case they read as 



|<|++(k,..) 



1 



1 



lc,|nvT"(k,c.) 



•-+0 >- 



|c,|nV'r"(k,c.) 



with 



70 
7 



A 

To 

(fcA)2(17-l)-7' 

To 

(fcA)2(f7-l)-27' 



l3No\c,\, 
l3M\c,\. 



(136) 
(137a) 

(137b) 
(137c) 
(137d) 



(137e) 
(137f) 



For kX <ti 1 and |r2fcA| <ti 1 the regular polarization function (the bubble) in the polar case can be approximated as 



(138a) 
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while in the ferromagnetic case they read as 

r('')++/i ^ 
^ (k,w) 



Cn,sn_|__|_ 



(Ml 



1 i n-i 

-- In ■ 



tt('')00/i ^ 



{^,lJ^kX2 n + r 

(AO 1 ^_^^^n-l + 2^4 



k\ 



27 



(138b) 

(138c) 

(138d) 



f7 + 1 + 2i- 



These results can be obtained with the Mittag-Leffler expansion of their spectral functions similarly as for the scalar 
case pH |. The approximation of the spin wave polarization functions can be obtained in a similar way, but due to 
their different symmetry properties compared with density polarization functions another limit is to be taken in the 
spectral function. The calculation is outlined in the appendix. The results are 



|cs|(n 

with 



(r)0+ 
+0 



n, 



(r)-0 



C- - C4 



2Cn + C4 



n 



r 



Co 

C+ 
C- 
D+ 



2r(|)f (II7) 

(ef)^ 

r(|)f-(i|o) 

(e)2 

K r(|)F(f|27) 

K 4r(|)f (||o) 

JC' 4r(|)F(f|27) 



B\2 



37r 



{k\f + O {{kXf) , (139a) 

(139b) 
(139c) 
(139d) 
(139e) 
(139f) 



and 7 is given by (|l37j ). 

In the case of linear dispersion SI ^ 1 can be assumed (while still satisfying the \^kX\ <^ 1 condition). For the 
ferromagnetic state we further restrict ourselves to ^J^ <^ \QkX\. We will see that this later condition is equivalent 
with c„/|cJ ^ 1, which is fulfilled in the applications. In this case the regular polarization functions with zero spin 
transfer (138) can be further approximated by 



|n(^)(k, 



i^ls? kXQ 



(140) 



both for the polar and ferromagnetic states. Since we are interested only in the low momentum behaviour of the 
dispersion curve Q will be searched as a power series of the (fcA) variable. 

Polar case, density mode (n ^ 0): The spectrum of these types of excitations can be determined by the poles of 
( |116a| ). The corresponding equation is 



1 - 3c„no - c„nc 



0. 



(141) 



Using Eqs. ( 137a ) and (14C) and substituting = a_i(fcA) ^ +0((fcA)°), the resulting equation for the a_i coefficient 
reads as: 



Solving this quadratic equation (and assuming, that the imaginary part is small) one arrives at 



n = ± 



(142) 



(143) 
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Returning to the more familiar variables one obtains 



OJ = ±\ Tr—k — I- TTTT^K 



M 



(144) 



for the beginning of the dispersion curve. 

Polar case, spin density mode (n = Q): The excitation spectrum is obtained as the poles of (116b). The corre- 
sponding equation is 



1 - 2csno = 0. 



With the use of Eq. (140) the resulting equation for Q. reads as: 

(Cf )2 nk 



l + 2i 



The solution is at 



or equivalently at 



. CsM 

^ = —1 — rrr^k 



(145) 



(146) 



(147) 



(148) 



for the beginning of the dispersion curve. Note that such a mode did not appear in the Bogoliubov approximation, 
because there Dzz is identically zero. 

Polar case, spin wave mode (n = ±1 ): The poles of the autocorrelation function ( 116cj ) determine the spectrum of 
these types of excitations. The corresponding equation is 



1 - 2c,no - CsHs = 



(149) 



Using Eqs. ( 137a ) and (|l40| ) and substituting 17 = a_i(fcA) ^ +0{{k\)'^), the resulting equation for the a_i coefficient 
in this case reads as 



all + 2i 



Whence 



or equivalently 



n = ±-\T,-i- ^ 



(150) 



(151) 



, INoCs. . CsM 

W = ±W — — K — I- — TTr-:k. 



M 



(152) 



Ferrom agnetic case, density, spin density mode (n — 0): The poles of the autocorrelation functions (122a) or ( 1221: ) 
or (122c) determine the spectrum of these excitations, which is given by the equation 



(1 - C„n„„) (1 - CsU^z) - CnCsU 



2 

nz 



0. 



(153) 



Substituting the approximations (1371:) and ( |140| ) and fl — a_i(fcA) ^ +C'((fcA)°), the resulting equation for the a_i 
coefficient reads as 



anCtsd'^i + ia{3as — 2an)a^__i — [as — an — &a^)a^i + iba ~ 0, 



(154) 



with a = /X, an,s — {S,n s/^)"^- the discussed temperature region as S> a„ S> a^. The solution of the equation 
can be calculated perturbatively, with results: 



las - a„ 1 .a 3a? + 2a^ 1 
il — ±\ -— — i- 



a^as k\ 2 a„as(as — a„) kX 



(155) 
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from which the frequency 



w ± 



M 



(156) 



We can see now, that the ■\/7 ^ |r2fcA| condition is equivalent with the 1 ^ 2c„/|cs| condition, which is wcU 
satisfied. 

Ferromagnetic case, spin wave mode with n = ±1: The excitation spectrum of the ferromagnetic spin waves is 
obtained from the poles of the autocorrelation function (122d). The equation to be solved is 



-I 



n: 



(r)04 



n: 



(r)-0 
0- 



n 



(s)04 



= 0. 



(157) 



Substituting the approximations Eq. (137c) and Eq. (139a) and that = oo + ©((fcA)^) leads to 
C- - C+ {k\f ' 



1 



To 
7 



7 



r 



(1 - ao)7o + 2Co + C+ + C- + ^(C_ - C+) + ^^—J^ 

7 7 



+ O{ik\f)=0 (158) 



The equation for 0{{kX)'^) is an identity, since 7 = 70 + f3\cs\{H^^ — H ). The ao constant can be determined 
from 0{{kXY) order terms 



ao = 



7(2Co + C+ + C_ + 70) -D++D- 
1' 



The frequency of the spin wave excitation then is a quadratic one 



2M* 



with the effective mass 



M* = — . 

ao 



(159) 



(160a) 



(160b) 



The Bogoliubov solution ( |106b| ) can be regained in the T —^ limit, since Cq, C+, C_ , D_ oc T, 7 and 70 tend to 

iVo(r = o)|c,|/fcBr. 

Ferromagnetic case, quadrupolar spin wave mode: Since the autocorrelation function (122e) is proper, so the 
spectrum of the quadrupolar spin waves is determined by the poles of the polarization function H^^. The singular 
polarization (119) has its pole at 



= 1 



27 

{kXf 



(161) 



The regular polarization function II^l ^ has a logarithmic singularity in the same place. The corresponding excita- 
tional frequency is 



2M\cs 



2M' 



(162) 



So me fin al remark is appropriate about the validity of the calculation of these last two modes. The approxima- 
tion ( 139a| ) is valid if (fcA)^/7 <^ 1. This is equivalent with the k <^ AM M.\cs\ / ft- condition, meaning that the 
approximation for the bubble breaks down near the transition temperature, where the magnetisation vanishes. 



VI. CONCLUSIONS 



In the high temperature phase density and spin fluctuations do not couple to each other and to the Green's func- 
tions. The latter appear only in the intermediate states in the perturbational expansion for the correlation functions 
describing density and spin fluctuations. Correspondingly one can experience independent excitation branches of 
one-particle and collective type. As discussed in the present paper in detail collective and one-particle excitations can 
hybridize in the Bose-Einstein condensed phase due to the symmetry breaking, which is different in the polar and 
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in the ferromagnetic phase resulting in different couphngs. In the polar phase, however such hybridization does not 
occur for spin modes characterized by spin transfers zero and ±2. The general results have been demonstrated in the 
RPA scheme. By treating the regular polarization contributions as perturbations, we obtained damping for a number 
of modes determined first in the Bogoliubov approximation |l^ , [l5[] . In the ferromagnetic phase for the transverse spin 
mode, whose energy has been shown to agree with the free particle kinetic energy in the Bogoliubov approximation, it 
is found that the eigenfrequency remains proportional to in RPA, but with an effective mass, which approaches the 
mass of the atoms in the zero temperature limit. Moreover the gap in the quadrupolar spin mode gets a temperature 
dependent correction, namely it is proportional to the magnetization which tends to the condensate density in the 
zero temperature limit and then the gap coincides with that of the Bogoliubov approximation. These transverse and 
quadrupolar spin modes are found free from damping. 

There are thermal excitations for which the bubble graphs can not be treated as perturbations as shown earlier in 
the scalar case [^,^. In the case of a spinor condensate one expects even a multitude of such excitations, already 
in RPA. Their investigation will be the subject of a forthcoming paper. 

The calculations in this paper has been made for a homogeneous system. In experiments with alkali atoms the gas 
sample is confined in an optical trap which can be modelled by a harmonic potential. The inhomogeneous nature 
of the trapped system results in that the algebraic equations presented here are to be changed to coupled integral 
equations but their main structure remains the same. We note also that experiments can be designed where the local 
speed of sound can be measured directly making the results obtained for a homogeneous system also experimentally 
relevant. 
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APPENDIX A: THE APPROXIMATIONS OF THE CONTRIBUTIONS OF THE BUBBLE GRAPHS 



In this appendix we briefly outline the analytical properties of the contribution of the bubble graphs encountered 
and outline the approximation (139a). Some of these results are known from earlier works |3l| . 

In the random phase approximation the self-energies of the internal lines are wavenumber and frequency indepen- 
dent, therefore the contribution of the bubble graph ( |10g| ), depicted in Fig. |l^ can be cast to the form 



n^y^r(k,zc.„) = - 



1 f d^q n°(ek+q + fiS'^'")-n"(£q-FftS'^''') 



with 



(27r)3 



(Al) 



AS 



(A2) 



The nondiagonal elements of Il^.'^^f'^ are zero, as stated below Eq. (10£), the same is true for the proper self-energies. 
In the appendix the automatic summation over repeated indices do not apply. 

Changing to retarded polarization functions with analytically continuing in frequency, the imaginary part plays the 
role of the spectral function: 



n 



'r)sr 



The imaginary part can be brought to the form 

^2 



{r)sr 



(k,c.') 



duj'. 



dz- 



2z 



dz- 



2z 



(A3) 



(A4) 
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with 







- 1 


2/1 




f 1 


7/ 






76 






A7 


= 7/ - 76 


= f 



A7 
A7 



and r2 defined in Eq. (136). Introducing the function 

2x 



R-y{z) 



dx, for 3mz > 0, 



which is to be continued to the whole complex plane, the function nrs''*'^(k, w) can be represented as 



(47r)2fi4fc^ 



yo 



Ry^ [z)dz - 



yi 



R^j {z)dz 



for Imuj > 0. 



(A5) 

(A6) 

(A7) 
(A8) 
(A9) 

(AlO) 
(All) 



The analytic structure of i?^(z) is important for the approximations, it can be found in Ref. |3l| ]. If the self-energies 
of the forward and backward propagating one-partic le lin es are equal, i.e. in thepolar case and for the n = mo de in 
the ferromagnetic case, A7 = and the expression ( All ) reduces that of Ref. |^l|. If the conditions above Eq. ( |l38| ) 
are fulfilled, the leading contribution of the Mittag-Leffler representation of R-y{z), reading as 



R^{Z) = 2V^(:{^) ~ iTTZ 



27r 



with 



47ri 
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«n +7 _ fen+7 
Z + an Z + bn 



+ (fln - K) 



ir„e 



ifin = arctan- 



2n7r 



(A12) 

(A13) 
(A14) 



can then be used to arrive at the approximating formulas Eqs. (138), see Ref. 

If the self-energies of the forward and backward propagating one-particle lines of the bubble graph are different, 
the upper limits of the integral (All) will diverge when k goes to zero, see Eqs. (A5) and ([A6|). If one is interested 
in the long wavelength dynamics the asymptotic series of R^{z) can then be used, namely 



R.y{z) 



2 r(n+|)f 

r,2 



2n 



(A15) 



where r( s) is the Gamma function and F{s\'y) is the Bose-Einstein integral (127) With the help of Eq. (All) 
and Eq. ( A15), the approximation (139a|) can be obtained. 
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